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A short review of both theoretical and experimental aspects of the ra- 
diative return method is presented. It is emphasised that the method gives 
not only possibility of the independent from the scan method measurement 
of the hadronic cross section, but also can provide information concerning 
details of the hadronic interactions. New developments in the PHOKHARA 
event generator are also reviewed. The 3 pion and kaon pair production is 
implemented within the version 5.0 of the program, together with contri- 
butions of the radiative <j> decays to the 2 pion final states. Missing NLO 
radiative corrections to the e + e~ — > /i + /i~7 process will be implemented 
in the forthcoming version of the generator. 

PACS numbers: 13.40.Ks,13.66.Bc 

1. Introduction 

Precise hadronic cross section measurements are crucial for predictions 
of the hadronic contributions to a a, the anomalous magnetic moment of the 
muon, and to the running of the electromagnetic coupling {(Hqed) from its 
value at low energy up to Mz (for recent reviews look [21 |Sj ) . When 
using the scan method one usually needs new experiments to be performed 
with not negligible costs, while the radiative return method proposed al- 
ready years ago |U allows for extracting the information on the hadronic 
cross section basing on the measurements at the existing meson factories, 
profiting from their huge luminosities and excellent detectors. As it was 
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shown in El the radiative return method is not only a powerful tool in 
the measurement of <r(e + e~ — ► hadrons), but allows for detailed studies of 
hadronic interactions. Due to a complicated experimental setup, the use of 
Monte Carlo (MC) event generators 01101111111211131111, which include 
various radiative corrections ^5] is indispensable for both signal and back- 
ground processes. Some more detailed analysis of that subject can be found 
also in [THllrTj. 

This paper is aimed as a short review of the results obtained by the radia- 
tive return method. It presents also a further potential of the method and 
shows new developments in the PHOKHARA event generator important 
for an extraction of the a(e + e~ — > hadrons) (and other physical quantities) 
from the measured cross section a{e + e~ — ► hadrons + photons). 

2. Hadronic cross section 

The extraction of the cross section a(e + e~ — > hadrons) from the mea- 
sured cross section a(e + e~ — ► hadrons + photons) relies on the factorisation 

da{e + e~ — > hadrons + wj) = H da(e + e~ — > hadrons) (1) 

valid at any order for photons emitted from initial leptons, where the 
function H contains QED radiative corrections. It is known analytically, 
if no cuts are imposed on photons, at next to leading order (NLO) and 
has to be provided in form of an event generator of the reaction e + e~ — * 
hadrons + photons for a realistic experimental setup. The emission of pho- 
tons from the final state hadrons has to be controlled as well [111 I17j . with 
an accuracy which allows for an error small enough not to spoil the accuracy 
of the a(e + e~ — * hadrons) extraction. 

Comparing the scan method and the radiative return method one has 
to say that they are in many aspects complementary. It is due to the fact 
that many experimental systematic errors are completely different in both 
cases. The radiative return method has though the advantage that most of 
them are the same for all the values of the invariant mass of the hadronic 
system for which the measurement is performed. That is not true for the 
scan method, where for each energy one has to perform a separate analysis 
of many systematical errors (an energy calibration etc.). It is also important 
that using the radiative return method one can use machines, which were 
built for other purposes (cp and B factories) and one has to 'invest' only in the 
experimental analysis. In many aspects both methods encounter however 
the same problems, which have to be solved. The already mentioned final 
state emission has to be studied carefully in both cases, the photon vacuum 
polarisation with its fast varying behaviour nearby resonances has to be 
taken into account, higher orders radiative corrections have to be properly 
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implemented in event generators etc. Even if the details are different for 
scan and the radiative return method the main features remain the same. 

The most important hadronic channel, from the point of the hadronic 
contributions to a,u, mainly 7r + 7r _ , is an example of that complementarity. 
Very accurate KLOE measurement |18| provided an important cross check 
of the CMD-2 data an d even if the agreement is not excellent it has 
allowed to conclude that the disagreement between e + e~ data and the r data 
concerning the pion form factor is not of the experimental origin. Further, 
mostly theoretical work, will be required to solve that puzzle and one will 
have to find new sources of the isospin violation effects, which finally will 
explain that disagreement. 

Already now many new valuable physical information was obtained 
by means of the radiative return method. The BaBar measurement of 
the a(e + e~ -> 2vr + 2vr-, 2K+2K~ , K + K~ti + tt-) [U| are the most accu- 
rate to date results, with the 2it + 2-k~ mode also extremely important for 
hadronic contributions to a„ and olqed- The BaBar measurement of the 
a(e + e~ — ► 7r + 7r~7r°) has shown that the cross section around the u>" 
resonance is actually much bigger as compared to an old measurement by 
DM2 collaboration |22| . Results coming from BaBar on narrow resonances 
|23l 1211 I2U| are also very competitive to the ones coming from the scan 
method (for a review see [H]). Not to mention the forthcoming results 
|25U26| . which show still growing potential for the radiative return method. 



3. Not only the hadronic cross section - looking inside the 

hadronic interactions 

The radiative return method originally proposed for the hadronic cross 
section measurements |lj|7j can be used to much more detailed studies of the 
hadronic interaction. The first investigations along these lines were done in 
[3] , where it was shown that it is feasible to measure separately the nucleon 
form factors in the time- like region at B-factories. That measurements are 
important for the understanding of the experimental situation in the space- 
like region (for a review look [27] ) , where two different type of measurements 
lead to different results for the ratio of the magnetic and electric proton form 
factors. 

The nucleon electromagnetic current is defined by the form factors as 
follows 

J M = - ie . u(q 2 ) (V(Q 2 )7„ " Ml) v(qx) , (2) 

with electric and magnetic form factors = +F^ , G^ = A-tF^ ■ 
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The statistics is not a problem for a measurement at B-factories with 
hundreds of fb _1 accumulated luminosity as the integrated cross section for 
the event selection corresponding to lower curve in Fig. (angular cuts 
close to BaBar angular acceptance) is about 59.3 fb for protons and 125 fb 
for neutrons in the final state. The separation of the electric and magnetic 
form factors is also possible for quite a big range of the nucleoli pair invariant 
mass. It is particularly easy if one performs analysis in the nucleon pair 
rest frame as shown in Fig. where the proton polar angle distribution is 
plotted both for a model which predict the ratio of form factors in agreement 
with measurements using the Rosenbluth method (G P M = n p G p E , triangles) 
and a model which predict the ratio of the form factors with agreement 
with the measurements using polarisation method (squares). For details 
concerning both methods see a review article [22j and references therein. 
For the 4 GeV 2 < Q 2 < 4.5 GeV 2 one expects about 2000 events per 
100 fb -1 accumulated luminosity and clearly two parameter (G^j,G E ) fit 
is possible and its accuracy will be limited mostly by systematic errors. 
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Fig. 1. (a) The differential, in Q 2 , cross section for the reaction e + e~ — * ppj for 
two different sets of event selections, (b) The differential, in cos 6 (6 - proton polar 
angle in the proton pair rest frame), cross section for the reaction e + e~ — * ppj for 
two theoretical models (see text for details). 



Already now BaBar collaboration has preliminary results for the proton 
electromagnetic form factor measurements |25j . which when completed will 
allow for extensive tests of the theoretical models. 

Another example 0, very specific for the radiative return method at 
DAPHNE energy, is the study of radiative eft decays at KLOE (for present 
status of the experimental situation see PS]). The <p decay (<fi ~^ fo(~^ 
7r + 7r~)7) contributing to the reaction e + e~ — > tt + tt~j might in principle 
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cause some problems in the pion form factor extraction from the <r(e + e~ — * 
7r+7r~7) measurement. However, as it was shown in jB], the detailed studies 
of the charge asymmetries allow not only to control that contribution, but 
also to distinguish between different models of the radiative <f> decays. That 
is clearly seen in Fig. where the charge asymmetry for an event selection 
enhancing the FSR contributions is presented. The differences between pre- 
dictions coming from different models of the radiative (f> decays and also from 
scalar QED (see [0] for details) are sizable in the region of large and small 
values of the pion pair invariant mass (Q 2 ) and a measurement can easily 
distinguish between them leading to tests with unprecedented accuracy. 



0.5 
0.4 
0.3 
: 0.2 
0.1 


-0.1 
-0.2 



7T+7T 7 



-0.3 



1 


™</> 


fO" 
f0" 


f0 KK 
10 str." 
~io str." 


1 

model — 9 — 

aj, = it 

a 4> = 2 

nofO -■ - ■ 


50° k 0k± 
-..50 • 9, 


< 13 
< 13C 


3° 



































v 

1 








i 



0.3 



0.4 



0.5 



0.6 0.7 0.8 



0.9 



Fig. 2. The pion charge asymmetry for different radiative <j> decay models 



4. New developments in the PHOKHARA event generator 

4.1. PHOKHARA 5.0: vr+vr-vr and KK final states 

In the newly released version of the PHOKHARA Monte Carlo generator 
(5.0) three new hadronic channels were added: 7r + 7r~7r°, K + K~ and KK. 
For the K + K~ both initial and final state photon(s) radiation was taken 
into account, while for the 7r + 7r _ 7r° and the KK only initial state photon(s) 
emission was considered. The kaon hadronic current was adopted from 28 , 
while a detailed analysis of all existing data on the e + e~ — > 7r + 7r _ 7r° cross 
section was performed in The constructed model allows not only for an 
excellent fit to the cross section and a good description of two pion invariant 
mass distributions (see Fig. |3J), but many three- meson couplings were ex- 
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Fig. 3. The two-pion invariant mass distributions: BaBar data [21] (points with 
error bars) and PHOKHARA5.0 predictions [2^1 (filled circles) . 



tracted separately from that fit making possible predictions of various decay 
rates and cross sections 

4-2. Next to leading order radiative corrections to the reaction 

e + e~ — ► /i+/i~7 

The reaction 

e + (p 1 )e'(p 2 ) -> v + (qi)v~(q2h(k) (3) 

may serve as a luminosity monitoring process for the radiative return method. 
If this method is used one measures the ratio 

o-(e + e~ -* hadrons + 7) 
a(e+e~ - + 7) ' W 

and for the extraction of the <r(e + e~ — > hadrons) from the data the the- 
oretical knowledge of both processes is needed. Due to a complicated ex- 
perimental setup that piece of information has to be provided in a form of 
event generators and the NLO radiative corrections to both processes are 
indispensable to provide accurate theoretical predictions. Already in a 
part of the NLO(FSR) radiative corrections was implemented and the miss- 
ing parts of the generator consist of diagrams shown schematically in Fig. |1J 
All details of the calculations and implementation in the Monte Carlo event 
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Fig. 4. NLO contributions to the <r(e + e — > 7) missing in the PHOKHARA 

5.0 code. 



generator PHOKHARA will be presented in separate publications |3U| while 
in this paper few details and the present status of this project is sketched. 
The contributions from diagrams in Fig. 0^ are implemented using helic- 
ity amplitude method in an analogous way as it was done for two photon 
ISR contributions 9 . The contributions from Fig. Et>, even if in principle 
are analogous to the ISR virtual corrections calculated in ^3], have to be 
treated in a different way as the generator has to work also for DAPHNE 
energy where the muon mass is not a small parameter. The expansions used 
for the ISR, which shorten the final result enormously, cannot be applied 
and even if the whole result is known further work is required to construct 
formulae for calculation of the radiative corrections fast enough for a Monte 
Carlo event generator. In both cases the helicity amplitudes can be writ- 
ten as a combination of 14 terms if the gauge invariance and the current 
conservation is used. The way one writes the result is not unique and we 
have chosen to write explicitly the part proportional to the Born amplitude 
plus suitably chosen symmetric and antisymmetric coefficients. The ISR 
corrections read then 
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Misr ~ u(q 2 )^^v{qi) ■ • v(j>x)S% u(p 2 ) , (5) 

i=i 

where the coefficients F{ contain loop corrections. The FSR corrections have 
analogous structure. 

The functions read 

S£ = W P+e* - fp+k) r , St; = \H p.e* - fp.k] r , 
S%=H \Pi e * P2k - P2(* Pik] , 5^ = k" [pie* p 2 k - p 2 e* pik] p^_ , 
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s? = fv% - H 



Pit M . P2^ fi 

— rP2 + — rPi 

pik p 2 k 



Pie n P2^ fi 

— r^2 rPi 

Pik P2k 



s i2 = e rP2 rPi 

Pik z p 2 k 

S i3 = \Pi e * P2k~P2(* Pik)pt , 5f 4 = [pie* p 2 k-p 2 e* Pik]p^ , 

where p± = p\ ± /?2 and e is the photon polarisation vector. 

Within that parameterisation all the coefficients, but the F\, are free 
from ultraviolet and infrared singularities, as the Si has exactly the spinor 
structure of the Born (ISR) amplitude. Moreover most of the coefficients 
vanish in the massless limit and they are numerically important only for the 
configurations with the photon collinear to one of the initial leptons. 

The diagrams in Fig. consist of box and pentabox diagrams. The 
pentabox tensor integrals (up to the third rank), which has appeared in 
the calculations, were reduced in D-dimension to standard box diagrams 
(tensor integrals up to rank two) with method equivalent to the one pre- 
sented in , even if in this particular case the reduction is simple due to 
the symmetry of the integrals. Further reduction is done using standard 
Passarino-Veltman reduction to scalar integrals. The reduction of two re- 
maining pentabox scalar integrals to the box scalar integrals has introduced 
simple denominators (kpi — kp 2 ) and (kqi — kq 2 ), which have zeros in the 
physical phase space. That problem will be solved by means of the expan- 
sion of the resulting around the mentioned zeros. Again here the remaining 
problem is the size of the result after the tensor reduction and further work 
is required to produce formulae, which can be used within a Monte Carlo 
program. 



5. Summary 

A short review of experimental results obtained by means of the radiative 
return method is presented, with an extensive discussion of the theoretical 
basis of the method. The new version of the PHOKHARA Monte Carlo 
event generator (PHOKHARA5.0) is presented. The status of the work on 
the NLO radiative corrections to the reaction e + e~ — > p + p~j, a luminosity 
monitoring process for the radiative return method, is also outlined. 
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